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FINITE ELEMENT ANALYSIS OF 
AEROACOUSTIC JET-FLAP FLOWS 
By 

A.J. Baker and P.D. Manhardt 
Computational Mechanics Consultants 
Knoxville, TN 


SUMMARY 


A computational analysis has been performed on the steady, turbulent 
aerodynamic flowfields associated with a basic jet-blown flap configuration. 

For regions devoid of flow separation, a parabolic approximation to the 
governing time-averaged Navier-Stokes equations is applied, which renders 
solution amenable to a downstream marching technique. Numerical results are 
presented for the flow on the symmetry plane of a rectangular slot-nozzle 
planar jet flap geometry, including detailed prediction of flowfield evolution 
within the secondary mixing region immediately downstream of the flap trailing 
edge. Using a two equation turbulence kinetic energy closure model, the 
numerical results predict rapid generation and decay of large spatial gradients 
in mean and correlated fluctuating velocity components within the immediate 
wake region. Modifications to the trailing edge turbulent flow structure, as 
induced by a simulated porous surface treatment of the flap, are evaluated 
using a hybrid turbulence closure model. A numerical analysis of the recircu- 
lating flow within a representative discrete slot in the surface is evaluated 
using a complete two-dimensional, time-averaged Navier-Stokes equation set. 

The parabolic analysis for a smooth flap is extended in an introductory manner 
to a finite span three-dimensional jet-flap flow. The results of the study 
are presented in this report. 


INTRODUCTION 

The use of directed jet flows is common in the design of aerodynamic lift 
systems. Examples include the leading edge slat-trailing edge flap configura- 
tions characteristic of current transport technology, as well as the lower 
and/or upper-surface blown flap geometries considered for STOL aircraft. In 
each instance, high momentum flow is directed generally tangential to an aero- 
dynamic surface. Such lift augmentation systems result in noise sources being 
generated by flow interaction with the lifting surface and equilibration with 
the free stream. With the progress made in noise reduction of propulsion 
system components, the noise floor associated with the next generation of pro- 
pulsive lift systems may well be constrained by the aeroacoustics of the 
fundamental jet-flap flowfield. 



Experimental testing of elementary configurations has been employed to 
characterize the aeroacoustic sources associated withthe basic jet*^f lap flow 
in an upper surface blowing (USB) orientation. Gruschka and Schrecker (ref. 1) 
evaluated a USB geometry comprised of a rectangular jet issuing over a planar 
flap with sharp trailing edge, and compared measured noise intensities with 
free jet results. Their data bear out the U® law for free jets, first noted 
by Lighthill (ref. 2, 3), and determined a 6th power law for the flap cases. 

The secondary flow mixing region, immediately downstream of the flap trailing 
edge contained a dominant noise source. Reshotko et al ., (ref. 4) tested a 
small USB model having a deflected circular jet issuing over a wing section, 
to determine the acoustic efficiency of the wing as a noise shield. The 
results showed increased effectiveness with increasing frequency similar to 
the results of Hayden (ref. 5). An extensive investigation of aerodynamic 
and acoustic phenomena of a slot nozzle and variable length straight flap 
was performed by Patterson et al. (ref. 6). They measured free field acoustic 
response, reverberation chamber acoustics, and utilized hot film anemometry 
and flow visualization techniques to correlate noise with flow perturbation 
phenomena. For the geometry tested, a maximum sound power level (SPL) 
occurred for a flap length of approximately 10 slot heights. An instability 
condition appeared for this case, as verified by flow visualization and hot 
film data. Longer flap lengths were determined to produce noise levels 
closer to the free jet measurements. Becker and Maus (ref. 7) report results 
of a comprehensive experimental project on the rectangular slot nozzle- 
planar flap geometry similar to reference 1. Using near and farfield micro- 
phone locations and a cross-correlation technique, they determined two extrema 
in acoustic source strength within the secondary mixing region, one located 
directly adjacent to the flap trailing edge. Detailed mean and fluctuating 
velocity correlation measurements indicate the turbulence structure in both 
mixing regions is highly anisotropic, and that sharp peaks in turbulence 
quantities occurred immediately downstream of the trailing edge. These 
rapidly dispersed as the flow proceded into the wake. 

These results generally indicate that a large portion of the fly-over 
noise associated with a USB equipped aircraft will be generated within the 
primary and secondary mixing regions. For the latter, the turbulent mixing 
flow and resultant acoustic source distribution is strongly dependent upon 
the boundary layer flow immediately preceeding the trailing edge (cf., ref. 8). 
Hayden, et al.,(ref. 5, 9) evaluated a variable impedance flap surface to 
reduce the noise intensity associated with a USB configuration. Penalties 
encountered in aerodynamic performance of the early systems were significantly 
reduced in more recent configurations, employing a cavity-backed porous mesh 
surface (ref. 10), while retaining the favorable broadband farfield noise 
reduction of 3 to 10 dB over a wide frequency range and for large turning 
angle. Many additional studies on powered and unpowered configuration noise 
measuremehts, as a function of flow parameters are reported (ref. 11-18) 
including standard configurations, and various aerodymamic components such 
as blown flaps, cavities, d-type surfaces, trailing edge interaction, and 
three-dimensional effects. 
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The results of these studies confirm the dominance of the jet flap flows 
as noise sources. A theoretical analysis would be directed at characterization 
of the basic mechanisms, and would require detailed information regarding the 
associated turbulent flowfield structure, particularly in the wake region down- 
stream of the flap trailing edge. Assuming the appropriateness of time-averag- 
ing, such data is potentially determinable by numerical solution of appropriate 
subsystems of the governing Navier-Stokes equations. For attached aerodynamic 
flows, the appropriate system is the boundary layer equation set, the numerical 
solution to which is routinely accomplished using mixing length theory for tur- 
bulence closure and any of several available solution algorithms. For free- 
mixing shear layer flows, as occur in the primary and downstream secondary mix- 
ing regions, the boundary layer set coupled with a turbulence kinetic energy 
closure model and algebraic length scale, is appropriate for a symmetric geo- 
metry. The wake flow within the immediate vicinity of the trailing edge is 
significantly more complex, and a complete analysis in the general case would 
require use of the full Navier-Stokes equation set. Such analysis could be 
extremely expensive, however, and simplifications have been proposed. For 
example, Melnik and Chow (ref. 19) employ a matched asymptotic analysis to 
characterize the trailing edge flowfield in a triple deck structure for lami- 
nar flows, with extension to turbulent flows for a symmetric geometry (ref. 

19, 20). Various forms of the boundary layer equation set have been employed 
as well for symmetric geometries (c.f., ref. 21). Numerical predictions for 
turbulent flows have been started somewhat downstream of the trailing edge, 
where the velocity minimum moderated (c.f., ref. 22). 

The present approach is to establish a parabolic approximation to the 
Navier-Stokes equation set by employing an order of magnitude analysis. The 
boundary layer equations are a simplified subset of the developed parabolic 
system, the use of which is not constrained to a symmetric geometry. A two- 
equation turbulence kinetic energy-dissipation function model is employed to 
close the developed system for turbulence phenomena. For aeroacoustic flows 
over flaps with sharp trailing edges, hence devoid of flowfield separation, 
the developed equation system can be marched directly off the flap surface into 
the trailing edge wake. Non-equilibrium turbulence phenomena within the imme- 
diate wake flow is allowed, such that local extrema in the turbulence phenomena 
can be predicted. The influence of a porous-acoustic treatment of the flap 
surface is simulated by appropriate boundary condition specification on the up- 
stream boundary layer flow. The influence in the resultant wake flow is then 
evaluated by direct numerical marching of the altered flow into the secondary 
mixing region. The validity of the porous surface simulation, regarding 
selected boundary condition equivalence, is evaluated by a complete Navier- 
Stokes numerical solution for the flow within the immediate slot vicinity. 

The developed parabolic concepts retain validity for non-separated three-dimen- 
sional flows over finite span planar flaps. 

The theoretical formulation of the aeroacoustic jet-flap flowfield model 
is presented. A brief overview of the basic flow illustrates how determined 
flowfield distributions may be employed in an aeroacoustic model. The required 
equation sets are presented including appropriate boundary condition specifi- 
cations as a function of the boundedness of the flow. A finite element algo- 
rithm is employed to cast the developed equations in form suitable for direct 
numerical solution. Results obtained using the COMOC computer program to solve 
these equations are presented to validate the developed concepts. 
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SYMBOLS 


sound speed; boundary condition coefficient 

Van Driest damping function 

wall porosity friction factor 

coefficient 

skin friction 

differential 

alternating tensor 

function of known argument 

drag force 

Froude number 

slot nozzle height 

boundary layer shape factor 

farfield acoustic intensity 

turbulence kinetic energy 

generalized diffusion coefficient 

differential operator; turbulence length scale 

differential operator; length 

finite element index 

Mach number; number of finite elements spanning R 

unit normal vector; nodes per element 

pressure; generalized parameter 

Stokes stress tensor 

generalized dependent variable 

generalized discretized dependent variable 


domain of elliptic differential operator 
Reynolds number 

acoustic source term; finite element assembly operator 
time 

Lighthill stress tensor 

velocity vector 

reference velocity 

friction velocity 

scale velocity 

observer distance 

Cartesian coordinate system 

Cartesian coordinate system 

friction velocity Reynolds number 

acoustic model parameter 

acoustic model parameter 

closure of solution domain R 

Kronecker delta; boundary layer thickness 

boundary layer displacement thickness 

turbulence dissipation function 

angle: boundary layer momentum thickness 

Karman coefficient (MLT) 

multiplier; turbulence sublayer constant (MLT) 
dynamic viscosity 
kinematic viscosity 



p 


density 

a-jj mean flow Stokes stress tensor 

T-jj Reynolds stress tensor; wall shear stress 

(f) finite element functional 

X generalized initial -value coordinate 

Tp streamf unction 

w turbulence damping factor; frequency; vorticity 

n solution domain 


Superscripts : 

e effective value 

T matrix transpose 

+ turbulence correlation function 

mass-weighted time-average 

time average 

'' unit vector 

^ mass-weighted fluctuating component; ordinary derivative 

* approximation 


Subscripts : 


CO 

e 


global reference condition 
freestream reference condition 
tensor indices 
jet reference condition 
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T s j » k j 

j 



non-tensor index 


m 

0 

t 

w 


finite element domain 
initial condition 
time derivative; turbulent 
wall reference condition 


Notation: 
{ } 

[] 

U 

n 


column matrix 
square matrix 
union 

intersection 
belongs to 



METHOD OF ANALYSIS 


Problem Description 

The general configuration for the aeroacoustic flowfield of interest is 
shown in Fig. 1, illustrating a source of high momentum fluid flowing over an 
aerodynamic surface subject to acoustic modification. The flow leaves the 
surface at a sharp trailing edge, tangent to the mean chord, in accordance 
with the Kutta condition, and proceeds to equilibration with the freestream. 
It is assumed the flow is essentially unidirectional and parallel to the Xj 
(curvilinear) coordinate as shown; hence Uj » U 2 , Uj where u.j is the 
velocity vector. 



Flap Deflection 
Angle ( 0 ) 

20 

60 

60 (Partially 
Detached) 


Figure 1. Schematic of Representative Jet Flow Patterns 
Over Wing/Flap Surface, Moo "" 0-8 (ref. 18) 


The point of departure for establishment of an acoustic model is the 
theory of Lighthill (ref. 2, 3). Based upon an exact analysis using first 
principles, Lighthill established that the partial differential equation 
governing propagation of sound in a homogeneous medium at rest is 

at2 0 ^ 9x.9Xj 


where a is the reference sound speed, and the solution domain is assumed 
0 

devoid of solid surfaces. Equation (1) is recognized as the wave equation; 
it possesses the retarded-time solution, expressed in terms of the perturbation 
to the mean density at the point x. , and the source strength distribution 
at y^ , in the form 


P 


32 


4iTa2 9x-3Xj 


^ij 


y,t 


X - y 1 

dy 

"0 J 

X - y 


( 2 ) 
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The noise source mechanisms are described by the Lighthill stress tensor, 
the right side of equation (1). They consist of the instantaneous 
convective accelerations and force terms as 


Tij = 


p"i“j 




- a^ p 6 . . 
0 1 j 


(3) 


The P-jj tensor contains the pressure and local viscous stresses and is 
expressed as 





y 




2 ^ 
3 3x,^ 



(4) 


A useful characterization of the terms in equation (3) is obtained by 
decomposition of the velocity field into mean time-averaged and fluctuating 
components as 


Hence, 


^i 


Ui + 


(5) 


Tij = p(u. + up(Uj. + uT) + 


^0 P j 


( 6 ) 


Evaluation of the second derivative of equation (6) is required; applying the 
continuity equation, the_^instantaneous source term of the acoustic equation 
in mean velocity field uj is 


9^T, 


S E 


8x .9x . 

* vJ 


^E T 


IJ’Jl 


o(u- .Q. . + Q. .Q. •) - 2u-d 
J’l I’J I’l i^»it 


u -U -0 • . 


+ 20- -(pul) . + 2(0. •pu''.) . + (pu'^ul) 


+ 


p. . . . 


a^ p . . 
0 


(7) 
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The density derivative terms are removed by a Galilean transformation. For 
small Mach number, temperature effects may also be neglected, which deletes 
the last two terms of equation (7). The source term that requires evaluation 
then becomes, in a moving reference frame 






+ 2(u. .pu") . + (pu"u") . . 


( 8 ) 


Previous analyses for elementary jet and wake flows have idealized the 
mixing layer by the assumption that only u^ is non-vanishing and that it is 
independent of Xj. This removes the first and third terms in equation (8) 
and simplifies the remaining summations. For the jet-flap flows of interest, 
however, the assumption that Ui= Ui(x 2 ) is inappropriate, especially at the 
flap terminus where large local accelerations can occur. The assumption on 
transverse mean velocity remains valid, however, which yields 


“ 20i_j(puj)_i + 2(Bi,iP>‘j),j * (pu:Uj), 




(9) 


The third term in equation (9) was originally analyzed by Proudman 
(ref. 23) using an isotropic turbulence model for free jets. The lead term 
was first identified by Mollo-Christensen & Marasimbo (ref. 24). In the 
terminology of Lilley (ref. 25) the third term is called the "self noise" 
due to its quadrupole nature. The first term is called the "shear noise" 
since the shear components are modified by the mean velocity derivative. 

The potential importance of the second term stems from the existence of the 
terminus of the trailing edge, a location experimentally verified to be a 
strong acoustic source. 


Substitution of equation (9) into (2), and utilizing a Galilean space- 
time transformation, yields the solution expressed in a reference frame moving 
with the flow as 


P - Po = 


^13 1 

= Ui •7^(pU‘T)+ = U- , .pu^ 

2iTa^M^ J 27ra2M J 


Lr"“o 


"j 


2ira^M^ 


u, i^(pup + — — ^(puCuT) 
1,1 3t J 4T7a^QM3 9t2 1 


dr 


(io) 


The farfield noise intensity has been determined from the variance of 
equation (10) in a fixed reference frame, as 


I(x) = — [p - Pq]. fp “ Pq 
PqL OJiL OJ; 


( 11 ) 


10 



The product is understood to include all possible tensor combinations. The 
overbar indicates time-averaged and equation (11) represents the noise 
intensity measured at observer location x due to all coherent sources. 

Evaluation of integrals in equation (11) is complex; modeling can be 
employed to express the covariances in terms of correlations of the turbulent 
flowfield. Equation (11) becomes a single evaluation in a uniform mean flow, 
with isotropic turbulence since only the self-noise term persists. Proudman 
(ref. 23) evaluated a simplified model using the concept of an eddy volume, 
beyond which significant coherence vanishes. He established the intensity at 
a point in the farfield for a moving reference frame as 

38p (k)^/2elx|‘^ 

I(x) = 5 

4ttM^6 (12) 

In equation (12), k is turbulence kinetic energy 


k = 



(13) 


which for isotropic turbulence is Uj Uj . The turbulence dissipation function 
e is defined as 



= 3v 


3Xk 9X^ 


(14) 


where v is the fluid kinematic viscosity. is an eddy convection factor. 

The derivation procedure of Proudman was applied to the mean shear noise 
term of equation (9) by Lilley (ref. 25). Under the assumption of isotropic 
turbulence, the farfield intensity is 


I(x) = C 


9u^" 


3Xr 


kJl5 


(15) 


where C is a constant and £ is a longitudinal turbulence length scale. 


Use of the concept for an axi symmetric free jet was proposed by Moon and 
Zelazny (ref. 26). For this case, the shear noise term, equation (9), is 
non-vanishing for j = 2, and several additional terms result from the' 
summation implied by the repeated subscripts in the self noise term. In the 
fashion of Lilley (ref. 25) and Csanady (ref. 27), noise was assumed radiated 
at two dominant frequencies. The time dependence in equation (11) was expressed 
in terms of these frequencies and appropriate eddy decay length scales. 
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Applying a directivity term for the free-jet shear layer flow, devised from 
geometrical acoustics theory (Csanady, ref. 27), and integrating over 
appropriately defined eddy volumes, the derived form for equation ( 11 ) was 

• ' 
j ^ ^se h h H 

87t2 |x| 2 nd - cos 9)2 

. 


9Uj 

“9r" 


^sh P" “sh 


^1 ^2 ^3 cos'^e + (u^up2 sin20 cos2e| 


2 fr 2 |x |2 [(1 - cos 6)2 + a 2 ^ ^/2 


dr 


:{16) 


The two terms radiate at the self noise and shear noise frequency 
respectively. The farfield intensity is an integral over the source field 
modeled in terms of turbulence parameters, i.e., components of the Reynolds 
stress tensor. For elementary two-dimensional or axisymmetric boundary layer 
and shear layer flows, the significant shear component of the Reynolds stress 
can be expressed as (cf., ref. 28, 29) 


■u^U2 = 


8 Uj 

3x^ 


(17) 


where C is a constant. The dissipation length scale is a function of 
k and c as 



(18) 


Substitution of equations (13), (17) and (18) then yields equation (16) an 
explicit function of the two-dimensional distributions of turbulence correla- 
tions and mean flow shear. 

The present focus is the more complex attached aerodynamic flow over and 
downstream of the terminus of a flap with a sharp trailing edge. Added 
complexity results from Uj becoming a function of both x^ and X 2 ; 

U 2 and U 3 remain negligibly small to first order. Several additional terms 
may assume importance in equation (9) as the j summations now range over 1 
and 2. Hence, both x, and Xj derivatives of Uj , as well as the turbulence 
correlations and several cross-product terms would result. A computational 
study of the basic flow geometry could initially focus on establishing detailed 
distributions of Uj and turbulence correlations, for example, k, e and 
and their derivatives. 
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The Flowfield Equations 

The basic aerodynamic character of the jet-flap flowfield has been 
illustrated. It is required to establish subsets of the governing Navier- 
Stokes equations that adequately describe the fundamental flow character and 
are also amenable to numerical solution. In Cartesian tensor notation, with 
summation implied for repeated latin subscripts, the non-dimensional form for 
mass and momentum conservation for flow of a compressible, single-species, 
isoenergetic perfect fluid is 


= It * “ 

• — 0 

LCpUj) = ^[p^j^i P'^ij “ = 0 

• • . J 


Cl9) 

C20) 


The dependent variables in equations (19)-(20) have their usual interpretation 
in fluid mechanics where p is mass density, Uj is the velocity vector, 
p is the static pressure, b is the body force. Re is the Reynold's number 
and Fr is the Froude number. The Stokes stress tensor, is defined in 
terms of the dynamic viscosity y as 




2n!!ii 

3 



( 21 ) 


The Navier-Stokes system, equations (19)-(21), becomes amenable to 
numerical solution techniques in a practical sense only after time-averaging. 
Employing the Reynold's decomposition (cf., ref. 30), define 


U, = U, + uj 


( 22 ) 


where u^ is the mass-weighted, time-averaged velocity 




p 


and u^ are the velocity fluctuations about the mean flow. 


to + T 


pu;: = 1 im i 

' J-yoo ' 


(pu^. - pu^. )dt = 0 


(23) 

By definition, 

(24) 
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L 


and 


pu,Uj = pa,Oj + puiui 


(25) 


The time-averaged equivalent of the Navi er-S tokes equations (19)-(20) becomes 


‘•(P) = » 

xJ 


(26) 


3(pu^.) 


L(pu,) = — ^ + # - (5.J. - puruj)] 

vl * 


= 0 


(27) 


where . is the time-averaged mean flow stress tensor. 


- P 
^ij " Re 


. 2 !!k 




6.-,- 


ax. 3 ax|^ ij 


(28) 


In eqn. 28, ia is the time-averaged dynamic viscosity, and the fourth term 
in the divergence is called the Reynolds stress tensor t.jj* 


T 


iO 


-pu.u • 
^ 1 j 


(29) 


We seek approximations to the steady-flow, time-averaged Navier-Stokes 
equations that yield adequate flowfield descriptions that are economically 
amenable to numerical solution using present day computers. One simpli cation 
is the parabolic approximation which can yield three-dimensional flow descrip- 
tions while requiring only two-dimensional computer storage. The three- 
dimensional parabolic Navier-Stokes equations (3DPNS) describe steady, confined 
or unbounded, viscous and turbulent flowfields wherein: 

1) a predominant flow direction is uniformly discernible. 


2) only in this direction are diffusion processes negligible 
compared to convection, and. 


3) no significant flowfield disturbances are propagated upstream 
against the predominant flow. 


Figure 2 illustrates the basic rectangular slot nozzle-planar jet flap 
configuration of interest, including representation of a finite element 
discretization and flap surface treatment, which is amenable to flowfield 
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, Fig. 2 Three-Dimensional Representation of A Rectangular 
Slot Injector - Planar Jet Flap Configuration 



characterization using the 3DPNS equation system. The predominant direction 
of flow is assumed pprallel to the Xj coordinate. The parabolic approxima- 
tion to equations (26)-(28) is accomplished by eliminating diffusion in this 
direction; hence, equation (28) becomes 


ij 


ii(l - iSii) 

3D. 

Re 

L®>‘J 


8u . 2 


6,.,- 


3x. 3x. 3 9x. IJ 

J I K 


(30) 


The mean flow unidirectionality assumption will also affect terms retained 
in the Reynold's stress model, as discussed in the next section. The sub- 
script bar notation denotes the index not eligible for summation, but is 
synonymous with the identical tensor index. 


The 3DPNS equation system contains, as a subset, the familiar two- 
dimensional boundary layer (2DBL) and two-dimensional parabolic Navier- 
Stokes (2DPNS) equations. Both these systems are employed to predict the 
jet-flap flow evolution on the symmetry plane of the three-dimensional 
geometry illustrated in Figure 2. The two-dimensional geometry is illustrated 
in Figure 3, including labeling of the primary and secondary mixing regions. 
For illustration, the 2DPNS equations in expanded form are 


L(p) 


3(pu,) 3(pu.) 

±_ + ^ 

3Xj 3X2 


3U] 


3u 


+ pu 


2 3x 


^ _ 1 3 


3 X, 


Re 3Xr 


3u, 


3x. 


3 

3x 


-pujuj 


_3 

3x 


2!- 


■pUjU2 


= 0 


3Uo r.- -1 r. 

L(pU2) = PU]L^'^ ^^2^"^ 3X2 " 


3u, 

i 

3x, 


3 

3Xj 



3 


■ -pu^U2 

3 X 2 

-PU 2 U 2 J 


0 


(31) 


(32) 


(33) 
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Fig. 3 Two-Dimensional Symmetry Plane Flowfield For Rectangular 
Slot Nozzle - Planar Jet Flap Configuration 



The dependent variables in the 3DPNS and 2DPNS system are the mean steady 
flow vector D-j and pressure p. For a two-dimensional boundary layer, the Xj 
pressure gradient is known from the freestream flow, and impressed across the 
boundary layer thickness. According to the same order of magnitude analysis, 
equation (33) vanishes identically to first order, yielding a second-order 
balancing of perturbation in pressure to the Reynolds stress normal component 
as 


L{p) = 



+ 


pu, 


)^2 


= 0 


(34) 


Closure of the 3DPNS system requires specification of • ; the turbulence 
model is discussed in the next section. 


The essential differential character of the 3DPNS momentum equation is 
initial -value in the x^ coordinate and bounday value on the (xg, X3) plane. 
Hence, given an initial distribution of u. , equations (27) -(28) are marched 
downstream parallel to Xj, and boundary conditions are imposed on the flap 
surface and at all lateral locations whereat the viscous jet flow merges with 
the assumed inviscid freestream. The starting solution plane and the location 
of boundary conditions specification is denoted in Figure 3 for the two- 
dimensional case. The boundary condition location for a three-dimensional 
case occurs everywhere along the outer extremity of the finite element grid 
illustrated in Figure 2. Correspondingly, on this closure segment, the jet 
velocity asymptotically matches the freestream value which is enforced as a 
gradient boundary condition. The velocity vector vanishes identically on the 
flap surface segment. Fig. 2-3, except if the flap is assumed porous, where- 
upon Uj takes on a specified non-zero value, i.e., U 2 (xj, 0, X3) = 

(Xj, X3). For the planar flap symmetry plane cases studied, the freestream 

pressure is uniform; therefore, to first order, the pressure is everywhere 
constant, and equation (34) provides a second-order estimate of pressure 
variation. The continuity equation (31) provides the freestream boundary 
condition for solution of equation (33) for transverse velocity in the 
secondary mixing region. For the boundary layer solution on the flap surface, 
equation (31) is solved directly for Uj as an initial -value problem in the 
X2 coordinate direction. Hence, a complete jet-flap flowfield solution 
requires a switching of equation solution procedure as the flow departs the 
flap trailing edge. 

The second simplification applied to the time-averaged steady-flow 
Navier-Stokes equations (26)-(28), for evaluation of a porous jet-flap flow, 
is reduction to two-dimensional space and transformation of dependent variables 
to a vorticity-streamf unction description (cf., ref. 31). In subsonic flows 
for which density may be assumed constant, equation (26) defines a divergence- 
free field, pu^- . From vector field theory, an equivalent expression on 
spaces spanned by Cartesian coordinates is 

(35) 

V.. 
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where e^.j|^ is the alternating tensor, J is the determinant of the metric, 

and is the streamfunction vector. Substitution of equation (35) into 
(26) yields an identity in zero. A useful transformation of equation (27) is 
accomplished by definition of the vorticity vector 

= 1. 

“i J jk (36) 

... ... . 

For two-dimensional problems, the sole non-vanishing scalar components for 
both and correspond to k = 3, in which case an elementary 2DNS 
differential equation set can be established. Denoting the Xs components of 
ijj and 0 ). as ijj and o), respectively, the compatability equation results 

K • K ' • ' 

frc)m substitution of equation (35) into (36) yielding. 


L(iP) 


9 1 9(f) 

9x . p 9X . 
J L 


+ 0 ) = 0 


(37) 


Taking the curl of equation '{27) to eliminate the pressure and substituting 
equations (35)-(36), yields the two-dimensional vorticity transport equation 
(cf . , ref. 31) , 


L(w) 



+ e 


3ki 9x, 


0) 


9ij) 

9x^. 



pufu;) 


= 0 


(38) 


where a . . 

ij 


is defined by equation (28). 


Equations (37)-(38) can be employed to evaluate two-dimensional transient 
low speed aeroacoustic flows wherein separation and recirculation are dominant 
features. Their use is appropriate, for example, for a detailed analysis of 
the flow departing a blunt trailing edge of a jet flap. For the present study, 
this equation set was solved for the recirculating flow within the immediate 
vicinity of a simulated porous slot on the jet flap surface, see Figure 3. 

The solution domain is shown in Figure 4, along with a representative finite 
element discretization. 
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Flow enters from the left, and the base of the slot region is assumed porous, 
hence a source or sink for mass. Equations (37)-(38) are both boundary value 
problem descriptions, hence boundary condition specification is appropriate 
about the entire closure illustrated in Figure 4. Any specified inlet/outlet 
velocity defines both i/; and to, using equations (35)-(36). Along the solid 
flap surface, ip is a constant and the no-slip boundary condition equivalent 
for to is (cf., ref. 31) 

0 ) = 

(39) 

where Xf, is the coordinate normal to the wall. Since equation (38) is 

also initial -value, the initial vorticity contour tOp is determined from a 
specified velocity distribution using equation (36). The initial distribution 
for rp is obtained from solution of equation (37) using tOp. 

The identified partial differential equations systems are potentially 
useful for determination of the turbulent aeroacoustic flowfields characteristic 
of the basic jet-flap geometry. It remains to establish a closure model for 
turbulence phenomena, to allow determination of the Reynolds stress tensor in 
terms of computational variables. 


Turbulence Closure Modeling 


The operation of time-averaging has introduced the Reynolds stress into 
the Navier-Stokes equations as well as the simplified subsystems identified 
for analysis of the jet flap flowfield. The primary requirement is for 
development of a closure model for the steady flow parabolic approximation, 
since the presented full Navier-Stokes analyses are restricted to flows 
dominated by wall damping. Using well known procedures (cf., ref. 32), the 
exact p artial differential equation description for the kinematic Reynolds 
stress -u^ u:[ in a steady mean flow is 

* w 


L(u'.uj) 
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^j"k 


9u. 
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9u^9u:; 
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(40) 
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9u:ju:: 

9x,, 




= 0 


where v is the kinematic viscosity, p/p. 
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Equation (40) is the departure point for development of a closure model. An 
additional differential equation for turbulence dissipation rate e is 
required; assuming the process is isotropic, 


„ _8ur3u: 

3^i3 ^ 


(41) 


The exact transport equation for dissipation function e is (cf., ref. 30). 


L(e) 
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3u^ 3u1 3u|; 
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k 3x^ 3x^ 
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(42) 


Equations (40)-(42) represent seven additional partial differential 
equations describing turbulence phenomena. However, this system is not 
closed since the third order correlations remain undefined. Additional 
differential equations could be established, but they in turn would involve 
undefined fourth order terms. Hence, modeling of third order correlations 
is invoked at a level of completeness, dependent upon the dimensionality and 
geometrical complexity of the physical system. For examp le. Launder et al 
(ref. 33) present closure in terms of all components of -u^uj. They document 

validity of the model for several cases including isotropic turbulence, free 
shear flows, elementary duct flows and flat plate boundary layer flows. 


In earlier work, Hanjalic' and Launder (r ef. 29). establish a closure 
applicable to thin shear flows where in only -UjUjis retained, and solved in 
combination with e and the turbulence kinetic energy k defined as 


k = iu'Tu" 
2 11 


(43) 


For the uni-directional, shear-dominated flows of primary interest, wherein 
Uj» u^, Ug, the contraction of equation (40) yields, after application of 

the parabolic approximation. 


L(k) 
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(44) 
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Equation (44) defines a new summation index convention appropriate for 3DPNS 
as, 1 < i > j < 3 and 2 < 5, <3. The corresponding form for the dissipation 
equation (42)~is 


L(e) 



'e 3x^. 


k e"^ 

• \J 


9e I 
3x. 

Jj 




e k 


3u, 

-1 1 




C" e" k- 
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1 _ 


= 0 


(45) 


In equations (44)-(45), the various constants are determined from approx- 
imate analyses and/or computer optimization (ref. 33). 

The next level of simplification involves specification of an effective 
turbulent diffusion coefficient v^. From first principles (cf., ref. 32), 
the effective diffusion coefficient must be of the form 


V. = CM Z 

(46) 

where C is a constant, V a scale velocity, and Z a scale length. For 
the turbulence kinetic energy-dissipation function two equation closure 
hypothesis (herein named TKE), V is taken as the square root of turbulence 
kinetic energy, equation (43). A dissipation length scale 5,^ defined in 
terms of k and e (ref. 34) is. 



The TKE closure hypothesis then specifies 


(47) 


V* 


= C.k^e"^ 


(48) 


Note that this is precisely the diffusion coefficient for turbulent kinetic 
energy, equation (44). Furthermore, upon summing the diffusion terms in 
equation (45), using equation (43) and assuming isotropy, equation (48) 
yields the diffusion coefficient for dissipation function as well. 
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To close the 3DPNS mean flow system, it is also required to model the 
correlation between the shear components of the Reynold's stress tensor and 
k, e and the mean velocity field u^- . Viewing equation (30), and neglecting 
dilitation, the required relation is assumed of the form 


-u-u„ 
1 Z 



£-1 



+ 



(49) 


The subscript bar indicates the index not eligible for summation. The elements 
of the correlation tensor are determined from a simplified analysis or 

experiment; the index range is restricted under the parabolic assumption. 


For the analysis reported herein, the 2DPNS and 3DPNS equation systems 
are closed assuming that in equation (49) is a diagonal Sensor. The 
overall effective diffusion coefficient can then be written as 


e _ 1 - 


y = 


Re 


y + pv. 


(50) 


The 3DPNS equation system for steady mean flow and turbulence closure, using 
the defined two-equation model and effective diffusion coefficient, then 
becomes 


L(p) = = 0 
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The tensor indices range 1 < i > j < 3 and 2 < a < 3 for 3DPNS. For 
symmetry plane analyses using 2DPNS, 1 < i,j < 2, 5, * 2 only. Hence, since 
i = 1 corresponds to the direction of predominant flow, diffusion is restricted 
to the plane transverse to the Xj coordinate, as required by the parabolic 
assumption. The recommended values for correlation coefficients for shear 
layer flows are given in Table 1 (cf . , ref. 33, 34). 


Table 1 

Coefficients in TKE Closure Model 


Variable 

Equation No. 

Coefficients 

^t 

(48) 

C^ = 0.09 

k 

(53) 

O 

II 

a. 

e 

(54) 

Prg = 1.3, C^ = 1.44, C^ = 1.92 


The boundary conditions for the mean flow equations have been described. 
Since the TKE equations {53)-(54) are also initial -boundary value descriptions, 
it is necessary to establish appropriate statements. Referring to Figure 3 
for example, the levels of k and e vanish in the non-turbulent freestream 
flow. Since the Reynolds stress hypothesis is valid only in regions where the 
turbulent Reynolds number is large, it is not economically feasible to enforce 
k and e to vanish at the flap surface. The alternative selected for these 
studies is to use boundary layer mixing length concepts to determine the 
distributions of k and e near the wall. Mixing length theory (MLT) 
expresses the correlation in equation (46) in terms of the predominant mean 
flow gradient and a length scale £ (ref. 35), 




8Ui 


9x, 


where £ is the mixing length 


£ 



0 < X2 < A6tc"^ 
X2 > X5k"^ 


(55) 


( 56 ) 


and 0 )' is the Van Driest function that accounts for the wall influence on 
velocity fluctuations. 
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OJS 1 - exp(-X2A“M 

In equation (56), Xg is the coordinate normal to the flap, 6 is the 
boundary layer thickness, and A and x are constants (0.09 and 0.435 
respectively). In equation (57), A is a function of many factors influencing 
flow phenomena near the surface including axial pressure gradient and normal 
mass flow addition. The form of Cebeci and Smith (ref. 32) serves to unify 
the many formulations as 
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(58) 


(59) 


All variables are time-averaged steady components, subscripts e and w refer 
to freestream and wall values respectively, A+ is a constant (25.3), and 
is the skin friction. Pressure gradient and mass addition affects are 
accounted for accordingly as 



"“le 

~3 

dx. 
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= u“^ V 
- w 


(60) 

(61) 


where Ug is the freestream axial velocity, is the specified transverse 
wall velocity, and u.j. is the shear velocity 


Ut - 


w 


IPJ 


The shear stress, Xyi is defined as 


8u, 


^w ^w^w 9x, 
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(62) 


(63) 
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The Ludwieg-Ti liman formula (ref. 36) yields 

_ 1 ~2 
2 ^le 


where Re^ is the Reynolds number based on boundary layer momentum thickness, 
and H h 6*0” where 6* is the displacement thickness (cf . , ref. 35). 

Equations (55)-(64) provide the formalisms necessary to determine k 
and e near a solid surface. These same concepts are employed to complete 
turbulence closure for the two-dimensional Navier-Stokes solutions for the 
pore slot recirculating flow, by identifying for equation (38) (cf., ref. 31), 
where ys is given by equations (50) and (55). 
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(65) 


Furthermore, through the dual definitions of turbulent effective viscosity, 
equations (48) and (55), and since the latter involves functions only of the 
mean axial velocity component Ui) which is either known or readily initial- 
ized, a means is established to initialize distributions of both k and e 
at the node points of a discretization. Since the developed TKE partial 
differential equations are initial -value, this information is required to 
start a solution. Additional comments on verification of this procedure are 
presented in the Appendix. 


The closure for turbulence phenomena is complete at the level of sophis- 
tication selected for these studies. The partial differential equations 
governing the flowfields of interest are now closed. All boundary conditions 
have been appropriately identified for partially and completely unbounded 
solution domains. The initial-valued character has been noted, and means 
established to initialize required distributions in terms of readily available 
data. Numerical solution of the developed system can provide the detailed 
distributions of mean flow and velocity fluctuation correlations required 
for a theoretical analysis. It remains to establish the numerical solution 
algorithm for these equations. 


FINITE ELEMENT SOLUTION ALGORITHM 


The desired form of the various partial differential equation systems 
governing the aeroacoustic jet-flap flows of interest are developed. Each 
is a special case of the general, second-order non-linear elliptic boundary 
value partial differential equation 


L(q) = 
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Herein, q is the generalized dependent variable, the tensor indices range 
2 < k,£ < 3 and 1 < i < 3, K is the diffusion coefficient, fj is a 
function~of its argument that specifically includes three-dimensional con- 
vection, p is a generalized solution parameter, and f 2 is the initial- 
value operator. The boundary condition statements for each of the dependent 
variables can be concisely expressed in the form 


S.(q) = 


ad), 





(67) 


i.e., the normal derivative of q is constrained by q and a parameter as 
determined by specification of the a(i)- An initial condition is required 
for q identified with each dependent variable as. 


q^x(O) (68 

The finite element solution algorithm is based upon the assumption that 
L(q) is uniformly parabolic within a bounded open domain n ; that is, the 
lead term in equation (66) is uniformly elliptic within its domain R, with 
closure 9R, where 


n = R’'[xo>x) (69) 

and xo ^ X- For the 3DPNS equations, x is associated with the Xj coordinate. 
For 2DPNS, it is time. Equation (67) expresses functional constraints on the 
closure of e br x[xo> x) > and the initial -condition specification, 

equation (68), lies on RU9R x xo' 

The concept of the finite element algorithm involves the assumption that 
each three-dimensional dependent variable is separable in the form 


q*(x.x^) = q^(x) q£ (X 2 .X 3 ) 


(70) 


The functional dependence in q 2 (x 2 s x^) is represented by a polynomial in Xj^. 
The expansion coefficients q^ can be most conveniently expressed in terms 
of the value of q*(x> Xjj^) at the nodes of the finite element discretization 

of R. Then, equation (70) takes the form 


q^(x»Xj2^) ^1 (^^ ^2 ^ ^3 
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where the polynomials 4>(x£) are known functions of Xj and Xg 
they are known, they can be differentiated analytically, e.g.. 
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(72) 


Hence, there is no need to establish difference formulae to approximate the 
differentiated terms in L(q). 


The finite element solution algorithm is established for the equation 
system (66) -(68) using the method of weighted residuals (MWR) formulated on 
a local basis. Since equation (66) is valid throught fi, it is valid within 
disjoint interior subdomains described by (x^. ,x) x [x ,x) > called 

finite elements, wherein UR^j, = R. The approximate solution f8r q within 

Rm x|xq,x)> called q* (x^. ,x) , is given in equation (71). Therein, the 
functionals are subsets of a function set that is complete on R^. 

The expansion coefficients Q|<(x) represent the unknown x~ dependent values 
of q* (x^. ,x) at specific locations interior to R and on the closure 3Rj^> 
called nodes of the finite element discretization of R. 


To establish the values taken by these expansion coefficients, require 
that the local error in the approximate solution to both the differential 
equation L(q^ ) and the boundary condition statement il(q*) for 3Rfrfl3R + 0, 
be rendered orthogonal to the space of the approximation functions. Employing 
an algebraic multiplier X, the resultant equation sets can be combined as 
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(73) 


where S is the mapping function from the finite element subspace R to the 
m m 

global domain R, commonly termed the assembly operator. The number of 
equations (73) prior to assembly is identical with the number of node points 
of the finite element R^. 


Equation (73) forms the basic operation of the finite element solution 
algorithm and of the COMOC computer program to be described. The lead term 
can be rearranged, and X determined by means of a Green-Gauss theorem: 
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For 8R^(18R nonvanishing in equation (74), the corresponding segment of the 
closed-surface integral will cancel the boundary condition contribution, 
equation (73) by identifying Aa*^) with K, equation (66). The contributions 
to the closed-surface integral, equation (74) whereat 3R^(19R = 0, can also 

be made to vanish. The globally assembled finite-element solution algorithm 
for the representative partial differential equation system then becomes 


S 

m 


ra ^ D 





(75) 


The rank of the global equation system (75) is identical with the total number 
of node points on RU9R for which the dependent variable requires solution. 
Equation (75) is a first-order, ordinary differential system for 3DPNS. For 
streamf unction in 2DNS, it is algebraic and the matrix structure is sparse and 
banded. Solution of the initial -value system is obtained by COMOC using a 
predictor-corrector finite-difference numerical integration algorithm. A 
banded Cholesky equation solver is employed to solve an algebraic equation. 

Solution is also required for the continuity equation (48) which is 
retained for boundary-layer flows. Since it exists in standard form as an 
ordinary differential equation, direct numerical integration yields the 
required solution at node points of the discretization. 
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COMOC COMPUTER PROGRAM 


The finite element solution algorithm is utilized, as observed in the 
previous section, to cast the original initial-valued, elliptic boundary- 
value problem description into large-order systems of purely initial-value 
and/or boundary value problems. The COMOC computer program system is being 
developed to transmit the rapid theoretical progress in finite element solution 
methodology into a viable numerical solution capability. COMOC integrates 
or equation-sol ves the discretized equivalent of the governing equation 
system. Initial distributions of all dependent variables may be appropriately 
specified or computed, and boundary constraints for each dependent variable 
can be specified on arbitrarily disjoint segments of the solution domain 
closure. The solutions for each dependent variable, and all computed para- 
meters, are established at node points lying on a specifiably nonregular 
computational lattice, formed by plane tri angulation of the elliptic portion 
of the solution domain i.e., RU8R. 

The COMOC system is built upon the macrostructure illustrated in 
Figure 5. The Main executive routine allocates core, using a variable 
dimensioning scheme, based upon the total degrees of freedom of the global 
problem statement. The size of the largest problem that can be solved is 
thus limited (only) by the available core of the computer in use. The 
precise mix between dependent variables and parameters, and fineness of the 
discretization, is user-specifiable and widely variable. The Input module 
serves its standard function for all arrays of dependent variables, para- 
meters, and geometric coordinates. The Discretization module forms the 
finite-element discretization of the elliptic solution domain and evaluates 
all required finite-element nonstandard matrices and standard-matrix mul- 
tipliers. The Initialization module computes the remaining initial para- 
metric data required to start the solution. The Integration module consti- 
tutes the primary execution sequence of problem solution, and utilizes a 
highly stable, predictor-corrector integration algorithm for the column 
vector of unknowns of the solution. Calls to auxiliary routines for para- 
meter evaluation (effective viscosity, Prandtl number, source terms, etc.) 
as specified functions of dependent and/or independent variables, as well as 
calls for equation solving algebraic systems, are governed by the Integra- 
tion module. The Output module is similarly addressed from the integration 
sequence and serves its standard function via a highly automated array- 
display algorithm. COMOC can execute distinct problems in sequence, and 
contains an automatic restart capability to continue solutions. A discussion 
on the functional design of COMOC is given in reference 37. 
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NUMERICAL RESULTS 

Numerical evaluation of the developed flowfield model can verify its 
utility by assessing accuracy of predictions. In concert with examination 
of a practical jet-flap flow, a test program was conducted to verify factors 
affecting accuracy of the finite element algorithm, as embodied in COMOC, for 
turbulent flow predictions. Both the MLT and TKE closure models were evalua- 
ted, and a discussion of results is presented in the Appendix. The following 
studies were conducted using discretizations and closure model combinations 
so identified to yield accurate results. 


Symmetry Plane Analysis of A Slot Nozzle-Jet Flap Flow 

The basic geometry involves interaction of a high momentum flow with a 
free-stream, over and downstream of the terminus of a planar jet flap, see 
Figure 2. Many experimental configurations have employed rectangular slot 
nozzles to form the jet flow, with aspect ratios (slot width to height) in 
the vicinity of 50:1. Important three-dimensional effects are then limited 
to the extremum boundary regions while the core flow approximates a two- 
dimensional character. Experimental data were taken (ref. 38) on the symmetry 
center-plane downstream of a slot nozzle-jet flap configuration of aspect 
ratio 60:1. This case was selected to evaluate predicted distributions of 
mean flow velocity and turbulence correlation. 

The basic experimental configuration and computational solution domains 
are illustrated in Figure 6. The jet flow is accelerated by the nozzle to a 
rsominal Uj = 120 m/s. Due to the associated favorable pressure gradient, 
the Uj profile at the starting plane of the solution is nearly uniform. 
Immediately downstream of the nozzle, the jet flow interacts with the free- 
stream within the primary mixing region, and a turbulent boundary layer flow 
develops adjacent to the flap which erodes the inviscid potential core at a 
rate different from the free shear layer mixing in the primary region. The 
flap terminates at a sharp trailing edge. Immediately thereafter, a secondary 
mixing region is engendered between^the jet boundary layer flow and the en- 
trained flow. The initially zero Ui on the flap surface is rapidly acceler- 
ated within the immediate downstream vicinity of the flap terminus. The large 
X 2 gradient of Uj associated with the turbulent boundary layer is conse- 
quently dissipated, and acts in the process as a strong source term for gene- 
ration of turbulence kinetic energy, see equation (44). Well downstream of 
the flap, the slot flow approaches a jet bounded by two free shear-layer mix- 
ing regions. 

The computational simulation of this composite jet-flap flow was accom- 
plished employing the 2DBL and 2DPNS solution options and the TKE closure 
model. The flowfield solution was initiated at the exit plane of the nozzle 
as noted in Figure 6. Since no experimental data were available to initialize 
the solution, a uniform profile for u.^ at the nozzle was assumed, and U 2 , 
k and e were started at zero levels. (This corresponds to the slug-start 
discussed in the Appendix.) The lower portion of the nozzle flow was assumed 
to develop on the flap as a boundary layer completely isolated from the primary 
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Fig. 6 Two-Dimensional Flowfield Specifications For Rectangular 
Slot Nozzle - Planar Jet Flap Configuration 


mixing region shear flow by the jet potential core. This was considered 
adequate, since experiment verified that the gradient-free potential core 
persisted well into the wake, and the far-field pressure was the constant 
associated with a free-jet flow over a planar surface. The free-shear layer 
mixing within the primary region was initialized using a step-profile in u^ 
as noted in Figure 6, where the upper Ui entrainment velocity was estimated 
from experimental data profiles at Xj/h = 6.6, one-step height downstream of 
the flap terminus (ref. 38). These boundary layer and primary mixing region 
solutions were marched downstream and then matched together at the trailing 
edge, x^/h =5.6. A new solution domain was specified to encompass the two 
flows, plus the lower entrained flow, and a 2DPNS solution initialized to pro- 
ceed into the wake region, see Figure 6. Again, the lower freestream Uj 
entrainment velocity was estimated from data, and was assumed laminar at the 
trailing edge. Boundary conditions for individual solutions within each 
domain are also noted in Figure 6. 

Using the 2DBL option, the flap boundary layer flow was started from the 
assumed uniform nozzle profile by numerical solution of equation (52) for Uj 
assuming laminar flow and the no-slip boundary condition, Ui(Xi, 0) = 0. 
Following the few integration steps required to establish derivatives, solution 
of equation (51) was initiated for computation of transverse velocity U 2 , 
assuming a non-porous surface. Following a few steps to allow equilibration, 
the developing laminar boundary layer was tripped turbulent by signalling 
computation of effective viscosity, equation (50), using the MLT model, equa- 
tion (55). The Uj and Q 2 profiles thereafter rapidly transform into the 
familiar turbulent profiles. The MLT solution was marched a short distance 
downstream, whereupon k and e initial profiles were computed using the 
dual definitions for vt equations (48) and (55). A restart of the entire 
solution was accomplished, and the turbulent boundary layer allowed to develop 

to the flap terminus using the TKE closure model, equations (53) -(54). The 

wall damping influence was retained within the TKE solution by over-riding 
the k and e levels, computed from the differential equation solutions, by 
those computed from MLT at all nodes lying inside the transitional layer. 

Illustrated in Figure 7 is the development of the turbulent boundary layer 

profile in terms of the shape factor H. For a fully developed, laminar in- 
compressible flat plate boundary layer, H - 2.6, while for a turbulent flow, 
1.3 < H < 1.6. The computed development spans the range. 

The two-dimensional shear layer computation within the primary mixing 
region was similarly initiated from a slug start. However, the computed Uj 
profile at Xi/ h = 5.6 exhibited a much larger potential core than did the 
experimental data of Xj/ h = 6.6. This indicates that the associated turbu- 
lent mixing within the blunt base region at the nozzle face was grossly under- 
estimated by the assumption of a thin shear layer. An accurate flow charac- 
terization in this region would require a complete Navier-Stokes solution, 
which could account for recirculation, as noted in the development. However, 
since the primary focus is on the secondary mixing region evolution, and since 
experiment shows the primary mixing region remains isolated downstream past 
Xi/ h = 6.6, the free-shear layer solution was simply continued downstream 
a distance sufficient to erode the span of the potential core to essential 
agreement with the data. The numerical profiles for u., u^, k and e at one 
slot-height upstream of this location were then employed to initialize the 
combined wake solution. 
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Fig. 7 Turbulent Boundary Layer Development Over Jet Flap 

The developed two-dimensional boundary layer and shear layer solutions 
initialized the 2DPNS solution on the combined domain at the flap trailing 
edge. The boundary condition specifications are denoted in Figure 6. Shown 
in Figure 8 are computed distributions of Uj at select downstream stations, 
shown in comparison to experimental data at x,/h = 6.6. Excellent agreement 
is illustrated, with the sole consequential difference at the wings of the 
Uj profile where the jet merges with the entrained flow. Shown in Fig. 9 
are corresponding computed distributions of turbulent kinetic energy k, equa- 
tion (13), in the primary and secondary mixing regions. Illustrated for com- 
parison is the measured normalized Xj component of the Re ynolds stress, 

-u" (ref. 38, Fig. 24). Assuming isotropy, k and Aij uj" would be 

directly comparable. Agreement is good within the secondary mixing region, 
wherein the initially small level at Xj/ h = 5.6 has been consequentially 
increased by the terminus of the flap. Considerably poorer agreement is 
noted within the primary mixing region, a direct consequence of the less 
accurate starting conditions as discussed. 

The illustrated agreement tends to confirm the validity of the wake flow 
initiation procedure as well as the appropriateness of the 2DPNS equation 
system. While the flow regions illustrated are important, the strong inter- 
action zone immediately downstream of the flap terminus is of primary interest. 
The extremum mean flow gradients exist therein, both in the Xi and Xj coor- 
dinate directions, and the corresponding generation rate of turbulence is 
extremal. Shown in Figure 10 are axial mean velocity u^ profiles at various 
stations downstream of the trailing edge. Note that the initial zero level 
on the flap is rapidly accelerated to produce the typical shear layer profiles. 
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Fig. 10 Mean Flow Velocity Profiles Downstream of Jet Flap Trailing Edge 


Since a mean flow velocity gradient in the direction continually persists* 
turbulence kinetic energy k is correspondingly generated well downstream of 
the local x. acceleration at the trailing edge. Corresponding computed dis- 
tributions of turbulence kinetic energy are presented in Figure 11 at identical 
stations downstream of the flap terminus. Following the initial maximum ac- 
celeration, the peak turbulence kinetic energy along a locus parallel to the 
flap in the wake continues to grow. The peak region broadens with distance 
downstream, and eventually generates a stepped peak. As a direct consequence, 
there results a pronounced overall increase in the level of turbulence within 
the flowfield due to the flap terminating. 

Numerical evaluation of spatial derivatives of Ui and k in the wake 
region could be employed in a noise model, for example, equation^(16) aug- 
mented for the more general case. The peak x^ derivative of Uj immediately 
downstream of the flap terminus is of the order 10’, and decreases rapidly 
to 10** one slot height downstream. The Uj derivatives in the x^direction, 
which contribute to the shear noise term, are also maximum at the trailing 
edge, rapidly decrease, and then continue to slowly^decrease as the flow pro- 
ceeds downstream. Calculated extremal values of 9ui/3X2 in the immediate 
wake are the order 10®, decrease to 10® at one slot height and .5 x 10® 
two slot heights downstream. The turbulent kinetic energy also contributes, 
and the extremal Xi derivative of k is of order 10’. Lateral x^ 
derivatives of k are also of order 10’, and persist well downstream of the 
trailing edge, as noted in Fig. 11. Recall that, as discussed, a noise model 
may utilize the eddy volume concept, the assumed-bounded region over which 
a non-zero correlation exists. A length scale for the eddy volume can be 
extracted from the computed turbulence parameters, see equation (47). Extremum 
dissipation lengths of .0075m and .0018m were calculated from the computed 
k and e distributions. They compare favorably with .0082m and ,0036m, 
as determined from experimental longitudinal and transverse space autocorrela- 
tion by Tam and Reddy (ref. 39). 


Acoustically Modified Planar Jet-Flap 

Treatment of an acoustically "hard" flap surface in the form of homoge- 
neous or discrete surface porosity is experimentally verified to alter far- 
field acoustic intensity (cf., ref. 5, 10). An aerodynamical ly-acceptable 
procedure is to replace the hard flap with a mechanically-formed mesh surface, 
through which mass flow can be induced by generation of modest pressure 
differences. Attenuation of far-field sound power level may primarily result 
from alteration of the local turbulent flowfield, in particular^the mean flow 
local shear stress distribution. Surface shear is basically 9ui/ 9xa, which 
is the dominant source mechanism in unidirectional shear flows (see source 
terms in equations (53)-(54)). It is also a mean flow contribution to the 
"shear-noise" source term correlation for the Lighthill equation solution, 
for example equation (16). The local value of shear at the flap trailing edge 
can be expected to significantly affect the rate of momentum defect attenua- 
tion within the secondary mixing region, immediately downstream of the trailing 
edge, see Fig. 6. Therefore, in this instance of interest, surface treatment 
appears to induce acoustic modifications by alteration of the local detailed 
flowfield structure. 
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Fig. 12 Influence of Porous Surface Treatment on Turbulent Flow 
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A computational simulation of the induced influence of surface porosity 
is particularly direct, since the equivalent boundary condition statement is 
a controlling mechanism of the numerical solution for the flowfield. While 
the actual mechanisms of pressure coupling may be rather complex, an elemen- 
tary extension of the two-dimensional configuration was evaluated to examine 
the fundamental phenomenon. The acoustic surface treatment can be assumed 
correlated by an elementary form of Bernoulli's equation. 


Ap = -cp 

^w w 


(76) 


In equation (76), Ap is the pressure difference across the porous surface, 
is the resultant induced colinear mass flux, and c is an experimental 
friction factor. The pressure difference exists between the exterior flow 
and a sub-surface cavity, which undoubtedly possesses a family of acoustic 
waves travelling at the characteristic cavity frequencies. Hence, one expects 
that Ap is a distributed function of x^. ; dependent upon the cited factors, 
the induced may be of variable sign and magnitude. 

The concept was evaluated using the detailed velocity field for the 
Bradshaw relaxing flow test case as discussed in the Appendix. This standard 
case was altered by specification of a discrete and cyclic distribution of 
Ap(xi), equation (76), as graphed in Figure 12a. The induced normal mass 
flux velocity, V^, acts directly as a boundary condition for solution of Uj, 
equation (51), and indirectly as a modification to the turbulence wall damp- 
ing function, equations (58)-(61). The flux period was approximately 0.02m, 
the wave-form a hat function with peak value = 0.001 (i.e., mass flow 

into the cavity:, and the Xi span of surface treatment was approximately 
0.15m. The flowfield variables were initialized from the standard case 
solution. Shown in Figure 12b is the computed Xj distribution of Uj , at 
the first finite element node above the flap surface (located at Xg/^ = 

0.0013 , where 6 is the local boundary layer thickness) in comparison to 
the standard case results. This modest efflux accelerates the local mean 
flow by up to about 10%, with a corresponding increase in BUj/ 9X2- The 
period appears equal to the applied pressure wave, and the phase lags by 
about one-quarter of the period. The effective turbulent viscosity, equation 
(50), computed at this node increased by approximately 8% at peak , which 

alters the corresponding value of k, equation (43), by about 15%. Reversal 
of the sign of Ap would induce deceleration of Uj by about the same mag- 
nitude, and the levels of v® and k would be correspondingly decreased. 
Elsewhere, away from the immediate vicinity of the surface, the computed Uj 
profiles were unaffected by the wall phenomena for the periodicity and 
amplitude evaluated. These results indicate that a very modest transverse 
mass flux, as induced by a pressure difference across a porous aerodynamic 
surface, can significantly alter the local detailed structure of a turbulent 
boundary layer flow. 


A corresponding influence on the local source mechanisms for acoustic 
phenomena could account in some part for measured alteration of farfield 


42 



intensity levels. In particular, the exact distribution of Ui and k at 
the flap terminus, as altered by transmission over the modified flap surface, 
could be important since this flow initializes the secondary mixing region. 

To assess influences, the discussed two-dimensional hard flap configuration 
of Schfecker and Maus (ref, 38) was computationally altered to induce a con- 
tinuous distribution of surface porosity. The wave form approximated a sine 
with amplitude Vw/Uoo = ±0.03 and period AXj/h = 0.04, where h is the 
slot nozzle height. In the downstream flap region, 5.0 < Xj/ h < 5.6, the 
standard case boundary layer flow was approaching fully turbulent with a shape 
factor H = 1.6, see Figure 7. This flow was rerun with the cited porosity 
distribution to determine the extremum induced modifications to H and the 
computed Ui directly above the flap surface. The results are summarized in 
Table 2. 


TABLE 2 

Porosity-Induced Jet-Flap Flowfield Modifications 


Case 

Velocity 

Shape Factor 
H 

No. 

Description 

Transverse 
V /u. 

w J 

Longitudinal 

u^/Uj 

1 

Standard (Hard) 

0 

0.092 

1.63 

2 

Influx (Soft) 

+0.03 

0.031 

1.87 

3 

Efflux (Soft) 

-0.03 

0.148 

1.57 


The period of the flow alterations agreed with the specified influence and 
lagged in phase as discussed for the check case. 

A 2DPNS solution in the trailing edge wake was completed to assess the 
influence of induced flowfield modifications upon evolution within the 
secondary mixing region. Solutions were initialized using mean flow velocity 
profiles for cases 2 and 3, that departed the furthest from the standard 
case 1. Table 3 summarizes the results in terms of longitudinal mean velocity 
and turbulence kinetic energy at two vertical (x^) levels and at several down.- 
stream stations. Even though the porous efflux case 3 is initialized with 
larger Uj , both the standard and influx cases produce higher u solution 
levels by Xj/h = 0.32. The standard case also produces extremal levels of 
k at this station, even though the initial levels for the efflux case 3 were 
five times larger. Note that by Xi/h 0.05 for the standard case, computed 
levels for both Uj and k are approximately 80% of those predicted at 
Xj/h = 0.32. Hence, the initial turbulence mixing phenomena within the 
secondary mixing region occurs directly adjacent to the flap trailing edge, 
and appears quite sensitive to the detailed structure of the turbulent boundary 
layer flow at the flap terminus. 
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TABLE 3 


Distribution of Mean Flow and Turbulence Velocities Within 
Initial Secondary Mixing Region as Function of Simulated Flap Surface Porosity 


Coordinates 

- x^-/h 

Mean 

Flow Velocity - 

“./Uj 

Turbulence 

Kinetic Energy -k/Uj^ 

Downstream 
(i = 1) 

Vertical 
(i = 2) 

Influx 
Case 2 

Standard 
Case 1 

Efflux 
Case 3 

Influx 
Case 2 

Standard 
Case 1 

Efflux 
Case 3 

0.0 

0.0014 

.031 

.094 

.148 

10-® 

10-® 

10'® 


0.0 

.0 

.0 

.0 

.0 

.0 

.0 

0.001 

0.0014 

.042 

.081 

.133 

.0006 

.0013 

.0063 


0.0 

.032 

.050 

.09 

.0005 

.0020 

.0044 

0.05 

0.0014 

- 

.344 

- 

- 

.0358 

- 


0.0 

- 

.320 

- 

- 

.0368 

- 

0.16 

0.0014 

.336 

.378 

.362 

.0300 

.0451 

.0376 


0.0 

.327 

.370 

.355 

.0300 

.0436 

.0374 

0.32 

0.0014 

.389 

.414 

.378 

.0331 

.0431 

.0382 


0.0 

.384 

.408 

.373 

.0335 

.0431 

.0382 










Recirculating Flow Within A Porous Slot 

The discussed calculations confirm that transverse mass flux through a 
porous flap surface can produce significant changes in the detailed structure 
of the attached turbulent boundary layer flow, the influence of porosity was 
exerted only indirectly on the Ui solution through the wall damping function 
and the Uj boundary condition. Consequently, a computational evaluation was 
made of the direct influence within the immediate vicinity of a porous slot 
by numerical solution of the complete two-dimensional time-averaged Navier- 
Stokes system, equations (37), (38), and (65). 

The test geometry corresponds to the flow over the first slot of the 
standard case as exhibited in Figure 12a. Both mass addition and deletion 
through the recessed base of the slot were evaluated. The solution domain 
geometry and boundary condition specifications are illustrated in Figure 13. 
Shown also is a representative discretization, illustrated with diagonals 
removed for clarity, consisting of 304 finite elements. The flow proceeds 
over the slot from the left, and upstream and downstream boundary conditions 
were established from the boundary layer solution using equations (35)-(36). 
Along the top of the domain, the gradient boundary condition on streamfunction 
allows the flow to respond to the cavity presence, i.e., the closure segment 
is not forced to be a streamline. Along the cavity base, was assumed a 

linear function of Xj with a maximum efflux/influx of = ± 0.015. 



Fig. 13. Finite Element Discretization And Boundary Conditions For 
Slot Region Recirculating Flow Analysis 
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The depth of the slot was approximately 0.5% of 6 the local boundary layer 
thickness. The flow Reynolds number based on 6 and Uoo is 10^. The slot 
Reynolds number based on shear velocity, equation (62), and slot depth, was 24, 
indicating that the contained flow was fully dominated by wall damping. Hence, 
the effective viscosity y® for all nodes within the slot was assumed laminar. 
MLT was employed to compute , equations (55)-(64), at all nodes above the 
original plate surface, as in the standard non-porous surface case. The top 
node row of the solution domain lies within the fully turbulent flow. 

Computed steady flow streamline distributions are shown in Figure 14a, b, 
for mass removal and mass addition through the slot base, respectively. For 
both cases, the far-field streamlines are computed concave downward indicating 
the effect of presence of the slot permeates the entire domain. Somewhat less 
concavity occurs for mass addition since identical farfield streamline levels 
are plotted for each case. The detailed flow structure in the immediate 
vicinity of the slot depends strongly on the sign of V^. For mass removal, 
Figure 14a, the boundary flow generally overshoots the slot and circles back 
along the base, establishing a closed circulation contour at the downstream 
extremity. Mass addition. Figure 14b, appears constrained to the slot region 
with emergence into the main boundary flow occurring at the downstream step 
face. Hence, within the assumptive constraints and for the specified boundary 
conditions, the more complete Navier-Stokes solutions further illustrate the 
character of flowfield alteration induced by a porous surface. These results 
are at best indicative however, since the flow through the cavity base was 
specified a priori, rather than being coupled to resonance phenomena within a 
sub-surface cavity. Nevertheless, they do confirm the potential capability 
to numerically establish detailed flowfield data of impact in an aeroacoustic 
analysis. 

An Elementary Three-Dimensional Evaluation 

The discussed numerical evaluations are constrained to flows on the 
symmetry plane of a three-dimensional flowfield. The three-dimensional flows 
associated with practical OTW configurations on airfoil surfaces, cf . , Figure 
1, are considerably more complex than amenable to analysis using the employed 
equation systems. However, for a rectangular slot-nozzle-planar jet flap 
geometry with quiescent freestream, such as that of ref. 38, an exploratory 
evaluation of the potential of a three-dimensional solution cam be established. 
Remaining within the boundary layer order of magnitude analysis, equation (34), 
and for the uniform freestream pressure field associated with rectilinear flow 
over a non-lifting surface, the mean flow between the jet symmetry plane and a 
lateral freestream, see Figure 2, can be approximated to first order as pre- 
dominantly boundary layer with negligible lateral velocity (us). Correspond- 
ingly, transverse mean velocity U 2 can be initially determined over the flap 
from the continuity equation solution equation (51), and thereafter by solution 
of the three-dimensional U 2 momentum equation, i.e., i = 2 in equation (52) 
and using equation (50). 
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a. Steady Flow Mass Removal Through Cavity Base 



b. Steady Flow Mass Addition Through Cavity Base 


Fig. 14. Computed Steady Flow Streamline Distributions For 
Turbulent Flow Over A Jet-Flap Slot. 
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A finite element discretization (less diagonals) of the left-half 
solution domain for the three-dimensional jet flap of Figure 2 is shown in 
Figure 15. Approximately 480 triangular elements are employed with non- 
uniformity specified to resolve the wall layer immediately adjacent to the 
upper and lower flap surfaces. The boundary condition constraints for the 
example solution are noted as well. As the flow departs the sharp trailing 
edge, the indicated wall boundary conditions are removed, as discussed for 
the two-dimensional solutions. A variant of the slug start was employed to 
initialize the solution field by interpolating the computed symmetry plane 
velocity distribution for u to a zero level at the lateral freestr'eam. 
Following initialization and equilibration of the computed Ug profiles, 
closure for turbulence was switched from MLT to TKE, and the solution marched 
downstream to a distance sufficient to smooth the solution field. The trail- 
ing edge was assumed to exist at this point and the solution restarted with 
the flap (i.e., boundary conditions) removed to simulate emergence into the 
wake. 


Shown in Fig. 16 is a surface representation of the computed equilibra- 
ted Uj distribution at the flap trailing edge. The grid imposed on the 
solution surface is identical to the employed discretization which serves to 
document appropriate refinement. The Qj velocity beneath the flap is assumed 
zero, and the additional grid detail therein has been omitted for clarity. 

The computed distribution of k at the trailing edge is shown in Fig, 17a. 

The centroidal spine is a consequence of the lateral derivatives of Uj , 
i.e., 8 Ui/3x 3 as produced by the appropriate contour in Fig. 16 which serves 

as a source contribution for turbulence kinetic energy, see equation (53), 

The peak directly above the flap is the three-dimensional equivalent of the 
results illustrated in Fig. 9. The extent and location of the flap is noted, 
and the lower discretization detail is again omitted. Fig. 17b illustrates 
the computed k distribution at AXj/h = 0.02 downstream of the flap 
terminus, as well as the lower discretization. The surface shape is unaltered 
everywhere except at the elevation of the flap terminus where a sharply spiked 
double peak has replaced the single peak illustrated in Fig. 17a. Hence, the 
existence of the trailing edge has resulted in rapid production of turbulence 
kinetic energy in agreement with the results of the symmetry plane evaluation. 
These three-dimensional predictions, while the result of a highly simplified 
analysis for an elementary geometry, do confirm the potential to extend the 
developed flowfield model to the three-dimensional configurations of practical 
interest. 
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CONCLUDING REMARKS 


Simplified forms of the Navier-Stokes equations for describing aero- 
acoustic flows over a basic jet-blown flap configuration have been established. 

A finite element formalism is employed to cast the identified ini tial -boundary 
value equation systems for steady, time-averaged turbulent flows into equiv- 
alent larger order systems of ordinary differential and/or algebraic equations. 
Numerical solutions were established using the COMOC computer program. 

Computed evolution of the turbulent flow on the symmetry centerplane of a 
rectangular slot-nozzle-planar jet flap geometry, and downstream of the sharp- 
edged terminus of the flap, compared favorably with experimental .data. The 
influence of a porous treatment of the flap surface on the detailed flow 
structure on and downstream of the flap terminus, was evaluated using two 
equation systems. An elementary extension to a three-dimensional flow con- 
figuration was evaluated. 

These results generally confirm the validity of the suggested approach 
to characterization of the turbulent aeroacoustic flows associated with 
directed-jet lift augmentation systems. In particular, use of the parabolized 
approximation to the full Navier-Stokes system appears appropriate for turbu- 
lent flows departing a sharp trailing edge of a planar flap. Extension to a 
curved flap surface requires development of a more comprehensive equation 
system, capable of computing pressure distributions in the plane transverse 
to the direction of predominant flow. This system should also be capable of 
predicting the entrainment induced by these lateral pressure gradients. 

Analysis of aeroacoustic flows over flaps with a blunt trailing edge using 
the parabolized equation systems is inappropriate in regions with flow separa- 
tion and recirculation. The present results indicate that the more complete 
analysis can be locally imbedded within a parabolic solution. In this instance, 
additional attention is required to adequately accomplish closure for turbu- 
lence phenomena in flow regions with small turbulence Reynolds number. 

Extension to these areas should render the developed concepts directly appli- 
cable to aeroacoustic flowfield determination for practical lifting configura- 
tions. 
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APPENDIX 


Initiation and Accuracy of Turbulent Flow Prediction 

A computational test program was completed to assess factors affecting 
solution accuracy of turbulent boundary layer flows predicted using the 
finite element solution algorithm. The standard two-dimensional boundary 
layer equations (2DBL) are a sub-set of equations (51)-(54), with equation 
(52) for j = 2 discarded and £ = 2 only,el sewhere. Equations (53) -(54) 
are solved for the turbulent kinetic energy (TKE) closure, using equation (50) 
to evaluate effective viscosity. For mixing length theory (MLT), the turbu- 
lence kinematic viscosity in equation (50) is determined algebraically using 
equation (55). 

The three evaluations required to attest to solution accuracy relate to 
verification of mathematical order-of-accuracy , turbulent flow solution 
initiation from an assumed mean-flow velocity profile, and the hybrid closure 
model employing MLT concepts within the wall layer to provide boundary con- 
ditions for the TKE solution. Regarding the first item, confirmation of a 
formal order-of-accuracy is currently evaluable only for laminar flows. The 
selected test case is laminar, incompressible flow at zero external pressure 
gradient, the well known Blasius similarity solution (ref. 35). The funda- 
mental error norm for a finite element solution is the energy norm (cf., ref. 
40), defined for the linear equivalent of equation (66)-(67) as 


E(q,q) 


3q 9q 


dr - X 


a^^^q^dxi 


3R 


(A.l) 


The 2DBL system is not linear, but the non-1 ineari ty exists in the lower-order 
convection terms only for laminar flow which would not constitute a quadratic 
contribution to equation (A.l). Assuming the validity of equation (A.l) for 
the 2DBL system, and for use of simplex finite element functionals, equation 
(71), convergence of the numerical solution is theoretically (ref. 40), 


E(q 


-q,5>q-q*) < CL±max|u 


— 


2 

"m 


where C is a constant independent of 
finite element spanning R , and uV is 


*'m 

the 


1 ' (A. 2) 

, the measure of the smallest 

second Xj derivative of u^ 
and assumed continuous on R. Hence, under discretization refinement, con- 
vergence should be quadratic in the energy norm. Since the Galerkin criteria 
for the finite element algorithm, equation (73), renders the error orthogonal 
to the approximation base, the exponent of two in equation (A. 2) can be 
numerically confirmed by measuring the finite element solution energy, equation 
(A.l) 


E(q*,q*) = 


M 

I 

m=l 


Rm 


9q* 9'q;; 
^m m 


dx - X 


,( 1 ) 


q*^ dx 


3R naR 

m 


(A. 3) 
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under discretization refinement. In equation (A. 3), qj!|^ is the finite 

element approximation functional, equation (71), ultimately evaluated in 
terms of the computed nodal distribution of the mean longitudinal velocity 
Ui(x, .X' ). Shown in Fig. A.l is the computed finite element solution energy, 
equation (A. 3), evaluated at a specific Xj station for 10 < M < 80, where 
M is the nominal number of finite elements spanning R(x^). The slope is 
almost identically two, indicating that the expected convergence rate is 
achieved by the computational embodiment of the algorithm. 


The laminar flow results were obtained using a uniform finite element 
discretization. However, a mandatory key feature for turbulent flow compu- 
tations is use of highly non-uniform discretizations to obtain satisfactory 
computational efficiency in concert with acceptable solution accuracy. 
Following numerical tests, solution speed and accuracy were both enhanced 
using a finite element discretization, variable according to a geometric 
progression as 


^m+1 


= + S 


m+1 

1 P 

j=2 


j-2 


1 :< m :< M 


(A. 4) 


In equation (A. 4), 11^+2 extremum nodal coordinate of and ri 2 

is the coordinate of the first node of the discretization, typically zero. 
Furthermore, p is the progression ratio and s is a scale factor that 
allows imbedding a given number of finite elements on R. Shown in fig. A. 2 
are graphs of discretizations using equation (A. 4) for several s/6 and p. 
Curves A, C and D illustrate uniform discretizations; the header shows the 
corresponding number of finite elements M spanning 6 and 35 , where free- 
stream boundary conditions are applied. Curve B illustrates a modestly non- 
uniform grid, suitable for laminar flow predictions. Curve E is the finite 
element discretization determined by numerical experiment to yield good 
solution accuracy for turbulent flow predictions in concert with minimal com- 
puter time. The finite element at the wall spans 6 x 10"^, yielding excellent 
resolution, yet only 28 finite elements are needed to span R = 36. 


Basic solution accuracy for a 2DBL turbulent flow was evaluated by com- 
parison of a zero pressure gradient, flat plate computation to the experimental 
data of Wieghardt (cf., ref. 36, Vol . II, IDENT 1400). Solution initiation 
was accomplished assuming existence of a uniform u^ profile at the plate 
leading edge and the MLT closure model. Boundary conditions at the pflate 
surface are^ = U 2 = 0; at freestream, x/6 = 3 , the vanishing gradient was 
imposed, 3u^/ 9 X 2 = 0. The initial distribution f(5r U 2 was zero since 
8 ui/ 9xx= 0 upstream of the plate. This specification corresponds to a 
"slug-start," the method employed to initiate solutions in the absence of any 
preferable alternative. Since the turbulent Reynolds number of the Wieghardt 
case is initially small, transition from laminar to turbulent flow occurs 
over a finite span of x,. The results of computational experimentation with 
the switch from laminar to turbulent flow is shown in Fig. A. 3, in terms of 
skin friction obtained from wall shear by non-dimensionalization with ^Poo ’ 
see equation (63). The compa^^ soft experimental results were obtained from 
data using interpolation and the Ludwieg-Til Iman formula, equation (64). In 
each case, the flow was computed laminar to the selected transition point. 
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and no transition model was employed to alter the intermediate profiles. The 
family of computed results are bracketed by the data, and by = 0.6m the 
various methods are in essential agreement. Corresponding comparison of 
computed u^ velocity profiles to data indicate excellent agreement, see Fig. 
A. 4. These results confirm the slug start solution initiation procedure for 
the selected planar jet-flap evaluations. 

The third requirement, to attest solution accuracy using the hybrid TKE- 
MLT closure model, was evaluated using as a comparison basis the experimental 
data of Bradshaw (ref. 36, Vol . II, IDENT 2400). Shown in Fig. A. 5 are com- 
puted Uj velocity profile distributions, illustrating the agreement with 
data attainable using the MLT closure model. Overall, the comparison is quite 
good, although diffusional processes within the flat mid-range appear high as 
evidenced by the computed results uniformly exceeding the data. These differ- 
ences diminish further downstream, but there is a corresponding trend to 
underpredict the first knee in the curve. Shown in Fig. A. 6 is the same com- 
parison to data with results computed using the TKE closure model. In the 
wall dominated viscous sub-layer, MLT was employed to compute near-wall 
boundary values of k and e , for numerical solution of equations (53)-(54). 
As previously mentioned, the MLT evaluation also yields the initial distri- 
butions for k and e. A vanishing normal gradient for k and e was 
enforced at the freestream, and computed agreement with data is comparable 
to the complete MLT run. Detailed differences do exist however, as illustrated 
in Fig. A. 7, on the multiple comparison bases of boundary layer integral 
parameters. The TKE colution does tend to overpredict both displacement and 
momentum thickness in comparison to the MLT solution. Otherwise the results 
are quite comparable and confirm the basic concept of the hybrid closure model 
for initiation of a parabolic solution in the wake downstream of a sharp 
trailing edge. 
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